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We propose and demonstrate a new read-out technique for a superconducting qubit by dispersively 
coupling it to a Josephson parametric oscillator. We employ a tunable quarter-wavelength super¬ 
conducting resonator and modulate its resonant frequency at twice its value with an amplitude sur¬ 
passing the threshold for parametric instability. We map the qubit states onto two distinct states of 
classical parametric oscillation: one oscillating state, with 185±15 photons in the resonator, and one 
with zero oscillation amplitude. This high contrast obviates a following quantum-limited amplifier. 

We demonstrate proof-of-principle, single-shot readout performance, and present an error budget 
indicating that this method can surpass the fidelity threshold required for quantum computing. 
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The readout scheme for quantum bits of informa¬ 
tion (qubits) constitutes one essential component of a 
quantum-information processor [Tj. During the course of 
a quantum algorithm, qubit-state errors need to be cor¬ 
rected; in many implementations this is done by quan¬ 
tum error correction, where each operation is based on 
the outcomes of stabilizer measurements that indicate 
the qubit errors. The stabilizers must therefore be de¬ 
termined in a “single shot” - without averaging of the 
output signals of repeated measurements on identically 
prepared qubits - with fidelity exceeding approximately 

99%®. 

The commonly used measurement scheme for a super¬ 
conducting qubit coupled to a linear microwave resonator 
does not, by itself, offer single-shot measurement per¬ 
formance. The qubit imparts a state-dependent (“dis¬ 
persive”) frequency shift on the resonator, which can 
be determined by applying a probe signal and measur¬ 
ing the reflected or transmitted signal, although only for 
weak probing, rendering an inadequate signal-to-noise ra¬ 
tio (SNR) HE]. 

Researchers have addressed the problem of insufficient 
SNR in essentially two ways. One approach is to feed 
the weak output signal into a following, parametric linear 
amplifier that adds only the minimum amount of noise 
allowed by quantum mechanics EH23- Another approach 
is to insert a nonlinear element into the system and ap¬ 
ply a strong drive tone, such that the resonator enters 
a bistable regime, hence enhancing the detection con¬ 
trast [MS- 

In this paper we propose and demonstrate a simplified 
readout technique in which a superconducting qubit is 
directly integrated into a Josephson parametric oscilla- 
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tor (JPO). We map the qubit states onto the ground and 
excited states of the oscillator, and demonstrate proof- 
of-concept, single-shot readout performance (SNR > 1). 
We obtain 81.5% qubit-state discrimination for a read¬ 
out time t = 600 ns; however, from the error analysis 
we infer a read-out fidelity of 98.7 ± 1.2%, taking into 
account known and reparable errors due to qubit initial¬ 
isation and decoherence (17.2 ± 1.2%). A realistically 
achievable qubit-relaxation time, Ti = 50 p,s, and a Pur¬ 
cell bandpass filter would reduce these errors from 17.2 % 
to < 0.5%, as well as shorten the required read-out time 
to t < 100 ns. The remaining errors, which are due to 
switching events in the oscillator (1.2 ± 0.3%), can be 
eliminated by improving the data-aquisition protocol - 
see Discussion and Supplementary Information. These 
qubit and detection improvements would bring the read¬ 
out fidelity to ~ 99.5%. 

Our readout scheme relies on parametric pumping of 
a frequency-tunable resonator by modulation of its in¬ 
ductance. The pumping amplitude exceeds the thresh¬ 
old for parametric instability, the point above which the 
resonator oscillates spontaneously, even in the absence 
of an input probe signal. This instability threshold is 
controlled by the state of the qubit, whose ground and 
excited states correspond to the nonoscillating and oscil¬ 
lating states of the resonator, respectively. In our mea¬ 
surement, the oscillating state produces a steady-state 
resonator field corresponding to 185 ± 15 photons, whose 
output we can clearly distinguish from the nonoscillating 
state when followed by a commercial semiconductor am¬ 
plifier, eliminating the need for a quantum-limited ampli¬ 
fier. Conceptually, this method can yield arbitrarily large 
contrast due to the parametric instability, and moreover, 
only requires a pump but no input signal. 

This readout scheme is well aligned with scalable, 
multi-qubit implementations. Parametric oscillators can 
be readily frequency-multiplexed |15] and allow for a sim- 



2 


plified experimental setup (compared to conventional mi¬ 
crowave reflectometry) without a separate input port to 
the resonator or a following parametric amplifier, and 
consequently also without additional bulky microwave 
circulators that would normally route the input and 
parametric-pumping tones. It is also possible to manipu¬ 
late the qubit via the flux-pumping line only, which fur¬ 
ther reduces the number of cables and interconnects. 


RESULTS 


The Josephson parametric oscillator (JPO). Our 

device consists of a quarter-wavelength (A/4), super¬ 
conducting coplanar waveguide resonator, shorted to 
ground in one end via two parallel Josephson tunnel 
junctions (JJs) — see Fig. 3 a )- The JJs form a super¬ 
conducting quantum interference device (SQUID), which 
acts as a variable Josephson inductance, Lj($, /(</>)) = 


$ 0 / 27t cos(7rd>/$ 0 ) \//q - / 2 (0) 


, where 1$ is the criti¬ 


cal current and $0 is the flux quantum. This inductance 
can be controlled by the external magnetic flux through 
the SQUID loop, <h(f) = + <h ac (t), and by the super¬ 

conducting phase difference across the JJs, 0(t), via its 
current-phase relation, I(t) = 7o sin 0(f). 

Time-varying modulations of $ and <fi - “paramet¬ 
ric pumping” - affect the resonator dynamics, albeit in 
rather different ways; moreover, the Josephson induc¬ 
tance is indeed both parametric and nonlinear. We ex¬ 
plain these differences in the Discussion section below. 
The resonant frequency of the JPO is parametrically 
modulated via the magnetic flux, $(i), which can lead 
to frequency mixing as well as parametric effects such as 
noiseless amplification of a signal, frequency conversion, 
and instabilities [THMSOl- 

The state of the JPO has a rich dependence on several 
parameters, some of which was studied recently, both 
theoretically [2TJ (22] and experimentally [8) [18j [20]. The 
equation of motion for the intra-resonator electric field 
amplitude, A , can be written as 


iA + eA* +SA + a\A\ 2 A + iTA = ^/2ToB(t). (1) 


Here e is proportional to the externally applied pump 
amplitude, <f> ac , which modulates the resonant frequency 
parametrically at close to twice its value, w p « 2u> r 
(degenerate pumping), and 5 = u] p /2 — u> r is the 
resonator’s detuning from half of the pump frequency. 
The field amplitude, A, and its complex conjugate, A*, 
are slow variables in a frame rotating at ui p / 2, and |H| 2 
is the equivalent number of photons in the resonator. 
The Duffing parameter, a, associated with a cubic 
field nonlinearity, arises from the nonlinear Josephson 
inductance. The linear damping rate has two compo¬ 
nents, r = Po + T r , where Po/27t = 1.02 MHz is the 
external damping rate, associated with the photon decay 
through the coupling capacitor, and T r /2tt = 0.30 MHz 
is the internal loss rate. The equation’s right-hand side 



Figure 1. Experimental setup and readout mechanism, 
(a) Schematic of the cryogenic microwave reflectometry setup. 
The transmon qubit (red) is capacitively coupled to the copla¬ 
nar waveguide parametric resonator (blue). The input and 
output flows of photons are denoted |B|“ and \C\ 2 , respec¬ 
tively, whereas the number of photons in the resonator is de¬ 
noted |H| 2 . The output signal is acquired using heterodyne 
detection of the amplified microwave signal. The components 
drawn in lighter gray are those that are rendered unneces¬ 
sary by the JPO readout method, thereby offering a simpli¬ 
fied experimental setup (see text), (b) Parametric-oscillation 
regions for the qubit ground state |0) (solid blue line) and 
excited state |1) (dashed blue line), respectively. These blue 
lines represent the instability boundaries, t = e t i,, where the 
number of steady-state solutions to Eq. 0 changes — see Eq. 
( |10| ) in Methods. The two panels on the right are measured 
[I,Q]-quadrature voltage histograms of the device output for 
the pump bias point indicated by the circles, revealing two dif¬ 
ferent oscillator states: I. Outside of the region of parametric 
oscillations, the resonator is “quiet” (|H| 2 = 0). II. Within 
the region, the resonator has two oscillating states (|H| 2 > 0), 
with a phase difference of n radians — see further Fig. [5] 


represents the input probe signal, such that \B(t)\ 2 has 
units of photons per second. The output flow of photons 
per second, |C(i)| 2 , is given by C{t) = B(t) — i\/2T a A. 
For low pumping amplitude, below the parametric 
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External dc-flux, F/tt 


Figure 2. Combined frequency spectrum obtained from 
qubit spectroscopy of the transmon qubit (in red) and through 
standard reflectometry of the resonator (in blue). The solid 
red and grey lines are fits. The dashed grey line, at resonator 
flux bias F = 0.185-7T, indicates the bias point at which we 
later demonstrate the readout method. Inset: Vacuum-Rabi 
splitting around the flux-bias point where the transmon fre¬ 
quency crosses that of the resonator. The minimum frequency 
splitting yields a qubit-resonator coupling goi/2n = 46 MHz. 


instability threshold, e < e th , this device works as a 
phase-sensitive parametric amplifier (JPA) for an input 
B(t) at signal frequency w s = lo p /2 [71 [IBHT5| [2?]. Note, 
however, that we keep B(t) = 0 in the measurements 
reported here. For a pumping amplitude exceeding the 
threshold, e > e th , spontaneous parametric oscillations 
set in — see Fig. |ljb) and Eq. (10) in Methods. The 
resonator field builds up exponentially in time, even in 
the absence of an input probe signal until it becomes 
limited by the Duffing and pump-induced nonlinearities 
and reaches a steady state mm- 

We connected a transmon qubit capacitively to the 
resonator [M — see Fig. flTa). The state of the JPO 
(oscillating or nonoscillating) can then be controlled by 
the qubit-state-dependent, dispersive frequency shift, 
X, which the qubit exerts on the resonator :25j [26] . 
When the JPO is being pumped above the threshold for 
parametric oscillation, with amplitude e and frequency 
detuning 5, then a change of qubit state effectively pulls 
the resonator to a different value of the detuning, outside 
of the region of parametric oscillations — see Fig. 0b). 
We denote the qubit-state dependent detunings by 
Jl°> = S — x and = <5 + x■ The resulting mapping 
of the qubit state onto the average number of photons 
in the resonator provides us with a novel qubit-state 
readout mechanism, which we exploit in this work. 


Characterisation of qubit and JPO. The de¬ 
vice and cryogenic experimental setup are depicted in 
Fig. 0 a). The sample is thermally anchored to the 
mixing chamber of a dilution refrigerator with a base 
temperature of 10 mK. The parametric A/4 resonator 
(in blue) is capacitively coupled to the transmission 
line (C c = 11.9 fF), yielding an external quality factor 


Q oxt = u r /2T 0 = 2 555. A transmon qubit (in red) is 
also coupled near this end of the resonator. 

The resonator output signal is amplified using a 4-8 
GHz high-electron-mobility transistor amplifier, with a 
noise temperature T)v = 2.2 K, followed by two room- 
temperature amplifiers. We detect the outgoing signal 
using heterodyne mixing. The signal is first downcon- 
verted to a frequency (w aF — w LO ) /2tt = 187.5 MHz; then 
the [I,Q]-quadrature voltages are sampled at 250MSs _1 , 
before they are digitally downsampled at a rate of 
20MSs -1 . 

We first characterise the transmon spectroscopically 
- see Fig. [2] - from which we extract the Joseplrson 
and charging energies, Ej/2n = 9.82 GHz and Ec/2n 
= 453 MHz, respectively. From the vacuum-Rabi split¬ 
ting, we extract a qubit—resonator coupling rate goi/27r 
= 46 MHz — see inset in Fig. [2] 

Next, we fit the frequency tuning curve of the resonator 
(with the qubit in the |0) -state) to the relation 

J r °\F)=u J r(F)~g 2 01 /A(F), (2) 

where F = 7r4> dc /$o denotes the static flux bias, nor¬ 
malised to the magnetic flux quantum. The effective dis¬ 
persive shift due to the qubit is 

A(F) (a (F)-Eg)’ ^ 

which, in turn, depends on the qubit-resonator detun¬ 
ing, A(F) = u a {F') - u r (F), with F' = F/ 8.88 + 0.58 
representing the effective magnetic flux of the transmon. 
Moreover, the qubit and resonator frequency spectra are 
well approximated by EMU 

Ua(F') « y/8Ej |cos(F')| F c - E C , (4) 


uj r {F) 


WA/4 

l + 7 o/|cos(F)|’ 


(5) 


where u} X /a/2t: = 5.55 GHz is the bare reso¬ 
nant frequency (in absence of the SQUID), and 
7 o = Lj{F = 0 )/L, = 5.3 ± 0.1% is the inductive 
participation ratio between the SQUID (at zero flux) 
and the resonator. The solid grey and red lines in Fig. [2] 
are fits to Eqs. ([2]) and Q, respectively. 


Single-shot qubit readout. We now demonstrate 
our method for reading out the qubit with the JPO. 
We choose a static-flux bias point F = 0.185 7r for 
the resonator SQUID, corresponding to a resonant 
frequency wi°^/27r = 5.218 GHz and qubit transition 
frequency w a /27r = 4.885 GHz — see dashed grey line in 
Fig. n|a). Consequently, the qubit-resonator detuning 
is A/27r = —334 MHz, and the effective dispersive shift 
is 2x/2ir = —7.258 MHz. We measured a Purcell-limited 
qubit relaxation time, X) = 4.24 ± 0.21 /us, and Ramsey 
free-induction decay time T| = 1.66 ± 0.32 /xs — see 
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Figure 3. Qubit readout by the Josephson parametric os¬ 
cillator. (a) Pulse sequence: The qubit 7r-pulse (in red), with 
Gaussian edges and a plateau of duration tv = 52 ns, is fol¬ 
lowed by a short delay, Td = 20 ns, before the pump is turned 
on at time t = 0. (b) The solid blue and red traces show 
the inferred photon number, \A \ 2 , in the resonator, with and 
without a prior 7r-pulse on the qubit, respectively. Note that 
the resonator latches, once it has entered into the oscillat¬ 
ing state, and remains there even if the qubit relaxes. The 
traces are the result of 10 4 averages of the raw data; the in¬ 
set shows a single instance of the raw data on the same time 
axis as the main plot. Prior to the sampling window of width 
t s = 300 ns, a delay t t = 300 ns is added to avoid recording 
the transient oscillator response. The hatched region around 
the average photon number represents our uncertainty, origi¬ 
nating from the amplifier gain calibration — see Methods and 
Supplementary Fig. 3. 


Methods, Supplementary Fig. 5, and Table 3. 

To operate the parametric oscillator as a high-fidelity 
qubit-readout device, we must be able to map the states 
of the qubit onto different states of the oscillator, which 
we must then clearly distinguish. We encode the qubit 
ground state | 0 ) in the “quiet” state (the empty res¬ 
onator) and the excited state 11 ) in the “populated” state 
of the resonator. Figure 0a) shows the pulse sequence 
for qubit manipulation and readout, and Fig. [3j)b) shows 
the resulting output from the JPO, operated with the 
pump settings <5 ^/T = —5.34, e/T = 3.56. 

The populated oscillator in Fig. 0b) contains 185±15 
photons. We obtained this estimate from a comparison 
between the probe-amplitude dependence of the resonant 
frequency and the expected photon-number dependence 
of the Duffing shift — see Methods and Supplementary 
Fig. 3. This number of photons should be compared to 


\A\ 2 = 200±3 photons, which is the solution to Eq. |TJ 
in the steady state {A = 0). 

In order to achieve such clear qubit-state discrimi¬ 
nation as in Fig. ib), we needed to make a judicious 
choice of flux-bias point, F, to mitigate the effects of 
two nonlinear shifts of the resonant frequency [20]. The 
Duffing shift dominates when F — ► ±7r/2, whereas a 
pump-induced frequency shift dominates when F —> 0. 
These shifts can move the resonator away from the 
proper pump condition, thereby effectively restricting 
the output power — see Methods and Supplementary 
Fig. 2. 

Moreover, the qubit—resonator detuning should be 
in the dispersive regime (A < 701 X in which the qubit 
state controls the resonant frequency of the resonator. 
Yet it must yield a sufficiently large dispersive shift, 
X > r (Eq. §, to produce clearly distinguishable output 
levels, corresponding to the |0) and |1) states. For our 
chosen flux-bias point, we identify the optimal pump 
settings by mapping out the parametric oscillation 
region as a function of pump frequency and amplitude 
— see Fig. 0a). 

An interesting feature is present within the left half 
of Fig. B a ) (where the populated resonator encodes 
11)). Here, when the qubit is initially in the |1) state, 
the resonator latches into its oscillating state for as long 
as the pump is kept on, and does not transition into 
its quiet state when the qubit relaxes, as one might 
have expected. This latching is shown by the blue 
trace in Fig. 0b). We attribute it to the existence of a 
tri-stable oscillation state in EH, associated with red 
detuning of the above-threshold region for the | 0 ) state. 
When the qubit relaxes, there occurs an instantaneous 
shift of the pseudopotential for the amplitude A , from 
bi-stable (with two 7r-shifted, finite-amplitude states; see 
Fig. 0b), panel II.) to tri-stable (with one additional 
zero-amplitude state). The field’s initial condition at 
the time of this shift, A 7 ^ 0, causes the resonator 
to maintain its oscillating state. A separate study 
of this latching feature will be reported elsewhere. 

We evaluate the obtainable state discrimination by 
collecting quadrature-voltage histograms at every point 
within the two regions of parametric oscillations in 
the [5, e]-plane — see Fig. 0b). We choose the pump 
operation point /T = —5.34, e/T = 3.56, indicated 
by the black circle, and show the characterization in 
detail in Fig. 0 In this point, the state discrimination 
has reached a plateau around 81.5%. Each histogram in 
Fig. 0a-b) contains in-phase (Vj) and quadrature (Vq) 
voltage measurements from 10 5 readout cycles, with 
each measurement being the mean quadrature voltage 
within the sampling time t s (blue window in Fig. 0 . 
We project each of the 2D-histograms onto its real axis, 
and thus construct lD-histograms of the Vj component 
— see Fig. 0c). We can then extract a signal-to-noise 
ratio, SNR = |/U| 1 ) — M|o) I /(o'! 1 ) + cr|o)) = 3.39, where /i 
and er denote the mean value and standard deviation, 
respectively, of the Gaussians used to fit the histograms. 
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Figure 4. Parametric oscillations, and state discrimination, (a) Output field of the resonator when the qubit is in its 
ground state |0) (top panel) and excited state |1) (bottom panel), (b) Contour plot of the state discrimination within the two 
parametric oscillation regions. The black circle in the left region, located at <5^ /Y = —5.34, t/Y = 3.56, represents the bias 
point used throughout the rest of the analysis and in Fig. §b). The state discrimination in this point is 81.5%. 


The peak separation of the histograms gives a confidence 
level of 99.998% for the readout fidelity. The peak 
appearing in the center of the blue trace arises mainly 
from qubit relaxation prior to and during the readout. 
We analyze this and other contributions in the next 
section, as well as in Supplementary Note 1 and Fig. 4. 

To extract the measurement fidelity from the his¬ 
tograms, we plot the cumulative distribution function 
of each of the two traces in Fig. [5jc) , by summing up 
the histogram counts symmetrically from the center and 
outward, using a voltage threshold, V th . From these 
sums, we obtain the S-curves of the probability to find 
the qubit in its ground state as a function of the voltage 
threshold value — see Fig. id). We define the fidelity of 
the measurement as the maximum separation between 
the two S-curves. 


DISCUSSION 

To evaluate the fidelity of the readout itself, as com¬ 
pared to the fidelity loss associated with qubit errors, we 
now present an error budget. From the histograms in 
Fig. [5jc) , we can account for 81.5% of the population, 
thus missing 18.5%. To understand the remaining con¬ 
tributions, we run a Monte Carlo simulation of the qubit 
population, consisting of the same number of 10 5 readout 
cycles as in the measured histograms. The simulation re¬ 
sults are binned in the same way as the measurements, 
using the Gaussian fits as boundaries, and taking into 
account the following statistics: (i) qubit relaxation and 


preparation errors, (ii) thermal population of the qubit, 
(iii) spurious switching events by 7r radians of the oscil¬ 
lator phase during readout (yielding a reduced sampled 
voltage), and (iv) peak-separation error due to the lim¬ 
ited signal-to-noise ratio. 

We find that the main contribution to the loss of fi¬ 
delity is due to qubit relaxation prior to and during 
the readout. From the measured relaxation time, T) = 
4.24 ± 0.21 /is, we obtain a fidelity loss of 11.6 ± 0.5%. 
However, this error can be reduced substantially (to 
<0.5%) by introducing a Purcell bandpass filter [28H30] 
at the output of the JPO; since the qubit is detuned from 
the JPO, this decreases its relaxation into the 50-0 trans¬ 
mission line. Such a filter would allow us to increase the 
resonator damping rate, r 0 , substantially reducing the 
readout time without compromising T \. This is shown in 
Supplementary Note 2 and Table 4. Note, however, that 
an increased resonator damping rate yields an increased 
width of the parametric oscillation region: consequently, 
the qubit-resonator coupling, < 701 , and detuning, A, need 
to be chosen accordingly to render a sufficiently large dis¬ 
persive frequency shift. 

From the simulation, we further attribute 4.5 ± 0.3% 
to qubit preparation errors. Another 1.1±0.4% can be 
explained from thermal population of the qubit; the ef¬ 
fective qubit temperature is T q = 45 ± 3mK. By adding 
these fidelity loss contributions due to the qubit to the 
measured state discrimination, we can account for 81.5 % 
+ 11.6±0.5% + 4.5±0.3% + 1.1±0.4 % = 98.7±1.2%. 

There are also errors introduced by the parametric 
oscillator itself: Switchings between the 7r-shifted oscil¬ 
lating states reduce the overall measured voltage. We 
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Figure 5. Quadrature voltage histograms of the parametric 
oscillator output, collected after digital sampling. The pump 
bias point was </ 0 ^/T = —5.34, e/T = 3.56 — see Fig. 0b). In 
panel (a), the qubit was in its ground state; in (b), a 7r pulse 
was applied prior to the readout pulse, (c) lD-histograms 
of the in-phase voltage component, Vj, from the quadrature 
histograms in (a) and (b). The black and white solid lines are 
Gaussian fits, from which we extracted a signal-to-noise ratio 
of 3.39. (d) Cumulative distribution functions, corresponding 
to the |0) and |1) states, obtained by sweeping a threshold 
voltage, Vth, from the center of the two histograms (V7 = 0). 
The maximum separation between the two S-curves yields a 
state discrimination of 81.5%. 


performed a separate control measurement that yielded 
2.4 ± 0.5 % switching probability, which translates into 
a maximal fidelity loss of half of that, 1.2 ± 0.25%. 
The switching rate of the parametric oscillator depends 
on many parameters, including damping rates and bias 
points; this error can therefore, with careful engineering, 
be decreased even further. We could, however, eliminate 
the effect of phase-switching events by using a rectifying 
detection scheme, e.g., a diode or a field-programmable 
gate array (FPGA), tracking the absolute value of the 
output field instead of its amplitude. 

The last and smallest contribution to the fidelity loss 
is the peak separation error, which accounts for the in¬ 
trinsic overlap between the histograms. However, this 
contribution is < 0.002% for our SNR of 3.39, and can 
therefore be neglected. For details on the error budget 
analysis, see Supplementary Note 1 and Fig. 4. 

By combining the above-mentioned improvements (re¬ 
duced qubit relaxation rate, optimised qubit manipula¬ 
tions and cooling, enhanced resonator output coupling, 
and rectifying data acquisition), the read-out fidelity 


could realistically reach ps 99.5%, limited only by the 
qubit relaxation. 

Finally, we demonstrate that the relaxation time of 
our qubit is not measurably afflicted by the pump - see 
Methods and Supplementary Fig. 5. Our measurement 
scheme is, in principle, quantum nondemolition (QND), 
see Supplementary Note 3; however, a proper experimen¬ 
tal and theoretical assessment of the back-action is out¬ 
side the scope of this work. 

Table[l]puts our results in the context of previous work 
on parametric and nonlinear Josephson amplification and 
detection circuits. 

A flux-pumped, parametric phase-locked oscillator 
(PPLO) was used as a following amplifier, also enabling 
sensitive qubit readout |5] . In our work, the qubit was 
directly coupled to the JPO, which simplifies the experi¬ 
mental setup by reducing the number of microwave com¬ 
ponents needed. Also, with a pumping amplitude be¬ 
low the parametric instability threshold, the flux-pumped 
JPA has been used to read out one qubit [7], as well as 
multiple qubits coupled to the same bus resonator [29]. 

There is another way of operating our device: instead 
of pumping the flux at u p ~ 2 u r , we can apply an alter¬ 
nating pump current (e = 0, B(t) ^0), now at a frequency 
close to resonance, u p ~ u r , and thereby directly modu¬ 
late the phase difference, </>. Both methods can provide 
linear parametric gain upon reflection of a detuned signal 
(w s ^ Up /2 and u s ^ u p , respectively). The flux-pumped 
JPA has a very wide frequency separation between pump 
tone and signal, because u s ~ u r ~ u p /2, which is a 
practical advantage since it makes the resonator’s entire 
instantaneous bandwidth available for amplification with 
no need to suppress or filter out the pump tone. More¬ 
over, the A/4 resonator has no mode in the vicinity of u p 
that the pump might otherwise populate. 

We emphasize that there are indeed two different phys¬ 
ical mechanisms in play, since flux and current pump¬ 
ing address orthogonal variables in the sense that $ = 
(ipi — v?2) < E’o/27t and </> = (</?i + <^ 2 )/ 2 , where tp\ and 
ip 2 denote the gauge-invariant phase differences across 
the two parallel JJs. This distinction is also evident in 
Eq. 0. The parametric flux-pumping term, eA*, mod¬ 
ulates the resonant frequency; it couples the resonator 
field amplitude and its complex conjugate, which can 
provide quadrature squeezing of an input signal and en¬ 
ables phase-sensitive parametric amplification; and for 
stronger modulation there is a parametric instability 
threshold into the JPO regime - see Fig. [T](b) . 

Current pumping by an input B(t), on the other hand, 
corresponds to an external force which directly con¬ 
tributes to the intra-resonator field A and drives its non¬ 
linear term a|A| 2 . For zero detuning, u s = u p , this is the 
driven Duffing oscillator which has no gain (it offers no 
phase-sensitive amplification); for stronger driving there 
occurs a dynamical bifurcation but no internal instability 
or parametric oscillations. 

Current-pumping with a moderate amplitude is used 
for linear amplification with the JPA ED E2, which 
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Device 

€ 
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B p 

# modes 

Ref. 

JPO ( * ) 

> Cth 

0 

0 

1 

This work 

JPA 

^5 e th 

^0 

0 

1 

El 

JPA 

^ e th 

^0 

0 

multimode 

m 

PPLO 

> Cth 

7^0 


1 

m 

JPA 

0 


7^0 

1 

m 

JBA (t) 

0 

0 

7^0 

1 

m 

JPC 

0 

7^0 

* 0 

2 

m 


TABLE I. Overview of different modes of operation for the 
various Josephson amplification and detection schemes. The 
variables refer to Eq. (IT]) , where e denotes the flux-pumping 
amplitude (at u p « 2cajj, and B s and B p denote alternating- 
current signal and pump amplitudes, respectively (at uj p ~ 
uj r ). The two readout methods marked with an asterisk (*) 
have the qubit directly integrated with the detector, whereas 
the other devices are used as following amplifiers. 

enabled, e.g., the observation of quantum jumps in a 
qubit jjjj. Current modulation is also used in the latching 
detection scheme of the Josephson bifurcation amplifier 
(JBA) [TU, TTJ [Tm [53 [53]. There, a higher-amplitude 
input strongly drives the Duffing nonlinearity near its bi¬ 
furcation point; the two qubit states can then be mapped 
onto two different resonator output field amplitudes. The 
JBA was used for quantum non-demolition measurement 
of a qubit, and in a lumped-element resonator |12j . in 
which a qubit-state sensitive autoresonance was observed 
in response to a frequency-chirped current drive. Yet an¬ 
other method is to couple the qubit to a linear resonator, 
which inherits a cross-Kerr nonlinearity from the qubit; 
current pumping of the resonator can then yield a strong 
output signal that depends on the qubit state mm- 
In conclusion, we have introduced a single-shot read¬ 
out technique for superconducting qubits — the Joseph¬ 
son parametric oscillator (JPO) readout. We demon¬ 
strated proof-of-principle operation, obtaining a bare 
state discrimination of 81.5 %. After correcting for known 
and reparable errors, this translates into an inferred read¬ 
out fidelity of 98.7 ± 1.2%, which by implementing a 
rectifying detection scheme can be further increased by 
1.2 ± 0.3%. With foreseeable improvements and opti¬ 
mization, this device would be an attractive candidate 
for implementing multi-qubit readout in the context of 
scalable error correction schemes. This fidelity and the 
readout time are both amenable to optimization. 

Our system integrates a parametric readout mecha¬ 
nism into the resonator to which the qubit is coupled, 
substantially reducing the number of components needed 
to perform single-shot readout in a circuit quantum elec¬ 
trodynamics architecture. Advantages offered by this 
readout technique include the potential for multiplex¬ 
ing and scalability with no need for signal-probe inputs, 
additional microwave circulators, or separate parametric 
amplifiers. As opposed to other integrated readout de¬ 
vices, our pump frequency is far outside of the resonator 


band and can thus easily be spectrally separated from 
other transition frequencies in the system. 
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METHODS 

Device fabrication. We fabricated our device on 
sapphire, using niobium for the waveguides and the 
transmon paddles, and shadow-evaporated aluminum for 
the Josephson junctions. To reduce the surface rough¬ 
ness prior to processing, the 2” c-plane sapphire wafer 
was pre-annealed at 1100°C for 10 h in an atmosphere of 
N2:02, 4:1, ramping the temperature by 5°C/min. The 
annealed wafer was then sputtered with 80 nm of Nb 
in a near-UHV magnetron sputter. The first patterning 
of the sample consists of a photolithography step to 
define alignment marks and bond pads, deposited using 
electron beam evaporation of 3 nm Ti and 80 nm Au. 
Next, the resonator, the transmon islands, and the 
pump line were defined in the Nb layer using a standard 
e-beam lithography process at 100 keV, and etched using 
inductively coupled plasma reactive ion etching (ICP 
RIE) in NF 3 -gas. 

The Al/AlOx/Al Josephson junctions forming the 
SQUIDs, used to terminate the resonator and for 
connecting the transmon islands, were then defined in a 
second e-beam step. After exposure, the 2”-wafer was 
diced into separate chips, using the exposed e-beam 
resist as a protective resist. Prior to the first evaporation 
step, the surfaces of the Nb films where cleaned using 
in-situ Ar-ion milling inside of the Plassys evaporator. 
However, due to the substantially different regimes of 
critical currents, Iq, required for the Josephson junction 
of the transmons and the parametric resonator, two 
sequential evaporations and oxidations were performed 
within the same vacuum cycle by rotating a planetary 
aperture mounted inside the evaporator load-lock, 
effectively shielding one half of the sample at the time. 
Finally, a post-deposition ashing step was done to clean 
the surfaces from organic residues. 

Finding the parametric oscillation threshold. It 

is hard to experimentally find the parametric oscillation 
threshold with good precision, when only considering 
the parametric oscillation region, Fig. [4^a) , whose 
observed shape gets smeared by the amplified vacuum 
noise. In this section we present an alternative method 
using a weak probe signal: we probe the parametrically 
amplified response as we sweep the pump amplitude 
across the instability threshold. 

We apply a probe signal on resonance, w s = w|°\ 
while applying a detuned pump signal, such that 
(uip — 2uj s )/2t: = 100kHz. The signal then undergoes 
degenerate, phase-preserving parametric amplification 
(red trace in Supplementary Fig. 1), while the para¬ 
metric oscillations are cancelled out since we measure 
the average amplitude of the field. The parametric 
amplification has maximum gain just at the thresh¬ 
old. We plot the magnitude of the reflected signal 
as a function of the pump power (at the generator), 
yielding an oscillation threshold P th = —10.8 dBm, as 
indicated by the dashed red line. As a comparison, 


we measure the output power of parametric oscillation 
(PO), for u} p —2io r = 0 and B(t) = 0 - see the blue trace. 

Limits of the parametric oscillation amplitude. 

As briefly discussed in the main text, there are two non¬ 
linear effects that move the resonator away from its pump 
condition, by means of their associated frequency shifts 

ESI, 


Aoj = —a\A\ 2 — /3T(e/r) 2 . (6) 


The Duffing shift dominates near flux bias F = ±7t/ 2; 
the Duffing parameter is approximated as 


a(F) 


7T 2 q; A /4Zo / 7 q \ 3 

Rk \cos(F) ) 


Oq 


\ 3 

7o \ 
cos(E )) 


(7) 


where Z 0 = 50 fi is the resonator’s characteristic 
impedance and Rk = h/e 2 is the quantum resistance. 

The pump-induced frequency shift dominates near F = 
0; it is approximated as 


P(F) 


T cos 3 (F) cos 3 (F) 

w a /47o sin 2 (F) °sin 2 (E) 


( 8 ) 


The resonator’s frequency tuning vs. F, Eq. (5) in the 
main text, is shown in Supplementary Fig. 2(a), for the 
parameters of our device, and Eqs. |7]) and © are plotted 
in Supplementary Fig. 2(b). This figure illustrates that it 
is essential to bias the system far enough away from the 
limiting points, F = 0 or 7r/2, such that neither frequency 
shift pulls the resonator too far from its pump condition, 
thereby severely limiting the attainable output power. 

The steady state solution of Eq. (1) in the main text 
yields an analytic expression for the expected number of 
photons within the region of parametric oscillations, 



which, for our analyzed bias point, amounts to 200±3 
photons in the resonator. From this number, we ob¬ 
tain a Duffing shift —o|A| 2 /27r ss —5.4 ± 0.3MHz (for 
a/ 27t = 27 ± 1.5 kHz per photon) and a pump-induced 
frequency shift —/3r(e/F) 2 /27r ss —0.64 MHz (for (3 = 
(7.5 ± 0.1) x 10 -3 ). 

The parameter (3 has the effect of skewing the para¬ 
metric oscillation region, yielding an expression for the 
thresholds plotted in Fig. 1(b), 


e 

r 


a/2/3 \ 


l-2/^± Jl-4/3 (p+ S - 


( 10 ) 


Calibration of attenuation and gain via the 
Duffing non-linearity. In this section, we present how 
we calibrated the gain of the amplifier chain, using the 
photon-number-dependent frequency shift of the Duffing 
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oscillator, —a|A| 2 , which we recall from the previous sec¬ 
tion. The frequency of the resonator as a function of 
input probe power takes the following form, 


W r (^s) — Coy(O) 


2aT 0 io(.Ps-Att-30)/io 

hu) r ( 0 ) 


( 11 ) 


where w r (0) denotes the resonant frequency with zero 
photons in the resonator, To and T are the external and 
total loss rates, respectively, and a is the Duffing fre¬ 
quency shift per photon — recall Eq. Q. Using Eq. ©> 
we can fit the extracted resonant frequencies as a func¬ 
tion of input probe power at different flux bias points, 
F, with the attenuation, Att, as the only fitting param¬ 
eter (since a can be extracted separately by fitting u>\u 
and 7 o — recall Eq. (5)). This is shown in Supplemen¬ 
tary Fig. 3, where the data for five different flux bias 
points are fitted to attenuations presented in Supplemen¬ 
tary Table 1. From these values, we obtain an average 
attenuation, (Att) = 127.5 dB, which can be compared 
with the installed 120 dB, indicating that we have a ca¬ 
ble loss of 7.5 dB at the measurement frequency. 

Moreover, from the same measurement we can also ob¬ 
tain an estimate for the gain of the amplifier chain by 
assuming that all the signal gets reflected when it is far 
off resonance with the resonator i. e. reflection coefficient 
15*!! | 2 = 1. Then, the gain is obtained from the relation 


G = \S n \ 2 + Att. (12) 


For the five gain estimates presented in Supplementary 
Table 1, we obtain a gain of G = 81.0 ± 0.37 dB, at our 
given bias point. The error bars for this gain estimation 
has two origins: ±0.17 dB from the residual of the linear 
fit to the gain values presented in Supplementary Table 
1, and another ±0.2 dB from the gain drift over time, 
which can be compared with our 91 dB of installed 
amplification. 


Calibration of the resonator photon number. 

From the obtained calibration of the gain of our ampli¬ 
fier chain, G, we can now calculate the conversion factor 
between our measured power on the digitizer and the 
number of photons in the resonator, using the following 
relation, 


Quantum coherence and readout nondestruc¬ 
tiveness. To study how the parametric pump strength 
affects the qubit’s relaxation time, we here present 
coherence measurements for the transmon. First, we 
calibrate a qubit pulse duration corresponding to a 
7r-pulse, using a Rabi measurement, where the pulse 
duration time is swept, for a fixed pulse amplitude. 
From the fit in Supplementary Fig. 5(a), a pulse length 
of tv = 52 ns was obtained, and the Rabi decay time was 
T rabi = 2.53 ± 0.15 fis. The histograms corresponding to 
the first 0.5/rs are plotted in Supplementary Fig. 5(b), 
using the same projective technique as for the histograms 
in Fig. 5(c) in the main text. Finally, we perform a set 
of Tj measurements for different pump amplitudes e/F, 
and compare these with traditional reflection readout, 
where we apply a weak resonant probe signal, but no 
pump ( B(t ) 0, e = 0). The fits to the relaxation times 
suggest that our readout is not any more destructive to 
the quantum state of the transmon than the traditional 
readout technique is. We note, however, that our 
extracted relaxation time is limited by the Purcell effect, 
yielding T\ ss [2r 0 (goi/A) 2 ] _1 = 4.11 /is. Also see 
Supplementary Note 3. 


|A | 2 


_ Ps - Pn _ 

2(ro/27r)fio;j CI> 10 G / 10 ' 


(13) 


where P s and P n denote our signal and noise power 
levels, respectively. We demonstrate this for the bottom 
panel of Fig. 3, where the resonator is probed at a 
frequency aj[ 0 ' > /2 tt = 5.212 GHz. The external damping 
rate is F 0 /27r = 1.02 MHz, and we calculate the back¬ 
ground power level from the end of the trace (when 
the pump is off). From the obtained SNR, the number 
of added noise photons can be estimated accordingly, 
|A| 2 /SNR 2 = 16.1 ± 1.3. 
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